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§ ; Abstract 

q^. If the CKM matrix element V u d, that can be derived from superallowed nuclear 

decays, neutron decay and pion beta decay, is used for a precision test of the unitarity 
of the CKM matrix, the combination of the present world data seems to indicate 
a small violation of the unitarity condition for the first row. While an accurate 
^ ■ calculation of the radiative corrections of 0(a) (RC) is crucial in order to determine 

the value of V u d as precisely as possible, the theoretical analysis has been limited in 
the past by the rather crude estimate of the effect of the hadronic structure. Only 
the contribution due to the axial current depends on the hadronic environment. 
We develop a strategy to deal with the influence of the hadronic structure on the 
decay properties of the simplest hadron, the pion, and calculate the contribution of 
the axial vector current to the RC, using a light-front model for the pion. Its qq 
bound state structure is well described by two parameters, constituent quark mass 
and confinement scale, that have been fixed by a comparison with the data. We 
take into consideration three different groups of two-loop diagrams, and derive their 
light-front representations. We discuss the associated zero-mode problem and show 
that the respective light-front amplitudes are free of spurious contributions. There 
is only a small model dependent uncertainty of the final result for the RC for pion 
beta decay. 



1 Introduction 



For three fermion generations, unitarity of the Cabibbo-Kobayashi-Maskawa 
(CKM) matrix requires the sum of the squared moduli of the first three elements to 
be equal to one: 

V 2 = \V ud \ 2 + \V us \ 2 + \V ub \ 2 = l. (1.1) 

A test of this property is of crucial importance since a violation of unitarity would be 
evidence for new physics, and the use of such a result to constrain possible extensions 
of the Standard Model would require a precise value of V 2 and its uncertainty. 

The unitarity sum V 2 critically depends upon the precise value of the matrix 
element V u d for the decays u — > dev and d — > uev. These quark level transitions 
give rise to superallowed Fermi beta decays, the decay of the free neutron n — > pev 
and pion beta decay n + — > ix^ev. In each case the measured rate can be used to 
determine the value of V u d, after radiative corrections (RC) and the effect of the 
hadronic environment have been separated out. 

A general formula for the RC of order a to the transition rates has been given by 
Sirlin [l|]. The total decay rate 1/r can be separated into the uncorrected expression, 
denoted by 1/t , and an an overall factor as 



1/r = 1/to 1 + 5) 



S = 



g (E ) + 3£n^ + A g 



a 



3(Qu + Q d )£n^ + 2C 
Ma 



;i-2) 



where Q u and Qd are the quark charges of u and d quarks. The Sirlin function 
g(E,E ) has been defined in as a function of the electron or positron energy E 
and represents the RC to the electron or positron spectrum in allowed beta decay. 
In the total decay rate 1/r it is replaced by the averaged value g(E ); E is the 
end-point energy of the spectrum. 

The correction terms of 0(a) consist of three distinct parts. The first two terms 
g{E ) + 3£n(Mz/M p ) represent the contribution of the vector current and are in- 
dependent of hadron dynamics. The Z boson mass Mz is a consequence of short- 
distance effects while the proton mass M p cancels in the sum of the two terms. The 
third term A g is a small asymptotic QCD correction term; A g = —0.34 [p], |3|. Fi- 
nally there is a contribution IniMz /Ma) + 2C induced by the axial vector current, 
where the logarithm is again the result of short- distance effects, with Ma acting as 
an effective low-energy cutoff (presumably roughly equal to the a>i meson mass), and 
2C stands for the remaining low-energy part. 

The value of Ma is uncertain; Marciano and Sirlin suggested a range AOOMeV < 
Ma < 1600MeK, while Sirlin M proposed an even wider range 

M a i/2 <M A < 2M al , (1.3) 
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with the central value at the ai meson mass M a \ = 1.26GeV. The quantity 2C 
is model dependent and has been calculated in the Born approximation in Refs. 
|| H using nucleon electromagnetic and axial form factors. For pion beta decay 
C = in the Born approximation since the axial vector current does not couple to 
a pseudoscalar meson. The resulting values for C are 

„ _ „ f 0.885 (superallowed and neutron decays) , , 

- C Bom - | Q (pi()n beta decay y (1-4J 

For superallowed beta decays there are additional nuclear structure dependent con- 
tributions to C which have been proposed and discussed in Ref. 0. 

The uncorrected decay rate l/r , defined by Eq.(1.2), still incorporates Coulomb 
corrections and Z-dependent radiative corrections of 0(Za 2 ) and 0(Z 2 a 3 ) for su- 
perallowed nuclear decays, and depends on hadronic form factors, which encode the 
effect of the quark structure of the decaying hadron. 

Recently the current status of V u d has been reviewed by Towner and Hardy [0, 
based upon the current world data for the three decay modes indicated above. To 
date, nine superallowed + — > + transitions have been measured to ±0.1% precision 
or better, and the result for V u d obtained from the average ft value is 

iKdl = 0.9740 ±0.0005. (1.5) 
From this value of V u d the unitarity sum, Eq.(l.l), becomes 

V 2 = 0.9968 ± 0.0014, (1.6) 

where the PDG 98 || recommendations for V us and V u b have been used in Ref.0. 
The value for the CKM matrix element V us determined from an analysis of kaon 
and hyperon decays is \V US \ = 0.2196 ± 0.0023, while the value for V u b is \V u b\ = 
0.0032 ±0.0008 and does not affect the unitarity sum at its present level of accuracy. 

According to the analysis of Towner and Hardy, the error bar associated with 
the value of V u d is caused mainly by the uncertainty in the RC (±0.0004) due to 
the prescription (1.3) for the effective low-energy cutoff and the uncertainty in the 
nuclear isospin symmetry-breaking correction (±0.0003), while the average experi- 
mental uncertainty is quite small (±0.0001). 

The problems associated with a precise treatment of nuclear structure effects 
can be avoided if the beta decay of free hadrons is considered instead. A survey of 
world data on neutron decay observables has been presented in Ref. and it has 
been noted that the derivation of the value of V u d from n decay is limited largely by 
the uncertainty in the overall average value of A = Qa/ov- However, there is a new 
result for the beta asymmetry obtained by the PERKEO II collaboration |J which 
leads to the value |A| = 1.2735 ±0.0021. This single value, combined with the world 
average for the neutron lifetime, leads to the following value for V u d [0 : 

\V ud \ = 0.9714 ±0.0015. (1.7) 
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The unitarity sum is then 

V 2 = 0.9919 ±0.0030. (1.8) 

The error given in (1.7) is three times larger than the error in (1.5) and is dominated 
by the uncertainty in the measurement of the beta asymmetry but, as in the analysis 
of the superallowed decays, still contains the contribution of the uncertainty in the 
RC. 

The results for V u d and the unitarity sum V 2 given in Eqs.(1.5)-(1.8) are con- 
sistent with each other and seem to indicate a substantial violation of the unitarity 
condition (1.1) for three generations. Moreover, they support the conclusion reached 
in Ref. 0, that the treatment of the effect of the nuclear environment in superal- 
lowed nuclear decays is reliable, with only a small error, and that there is no evidence 
that the unitarity problem can be solved by improvements in the calculation of nu- 
clear structure effects. 

In order to obtain more information on the unitarity problem accurate measure- 
ments of the pion beta decay observables would be of great importance. Like the 
superallowed nuclear decays pion beta decay is a pure vector transition and the ma- 
trix element of the axial vector current, which complicates the analysis of neutron 
decay, does not contribute to the lowest order amplitude. In higher orders both the 
vector and the axial vector parts of the weak current contribute. The expression for 
the radiative corrections of 0(a) is given in Eq.(1.2). Moreover, since the decaying 
pion is free, the nuclear structure dependent corrections that complicate nuclear 
beta decay are absent. Based on the lifetime || 

r exp = (2.6033 ± 0.0005) x 10" 8 s (1.9) 



BR = (1.025 ± 0.034) x HT 5 , (1.10) 
the value of V u d was determined in Ref. |7j to be 

\Vud\ = 0.9670 ±0.0161 (1.11) 

and the unitarity sum 

V 2 = 0.9833 ±0.0311. (1.12) 

The price to pay for the advantage of a simple theoretical analysis of pion beta decay 
is a large error in V u d due to the considerable experimental difficulty in measuring the 
7r branching ratio with a precision comparable to the one obtained in superallowed 
beta decays. However, there is a proposal for an experiment at PSI [|TTJ with the 
aim of making a precise determination of the pion beta decay rate. In the first phase 
of the experiment it is intended to measure the branching ratio with an accuracy of 
0.5%. The proposed experimental method was designed to finally achieve an overall 
level of uncertainty in the range of 0.2 — 0.3%. 



and the branching ratio 
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The decay rate for pion beta decay including the RC of order a is given by Eq. 
(1.2), where an approximate expression for the uncorrected decay rate has been 
derived long ago by Kallen [|l!| : 



1/ro 



G 2 F \V ud \ 2 / A 



307T 3 

15 



2M+ 



9e 
~2~ 



4e 2 



1 + yO^ 



A 5 /(e,A), 



A 2 



7 (M + + M ) 2 



;i.i3) 



;i.i4) 



with e = m 2 /A 2 and A = M + — M , where M + and M are the masses of ir + and 7r°; 
Gf is the Fermi coupling constant. Equation (1.14) includes the leading correction 
in an expansion in powers of A 2 / (M + + M ) 2 [|I|. The effect of the quark structure 
has been neglected entirely, and in Sect. 2 we shall study the error made by this 
approximation. In particular, we shall investigate the effect of isospin violation due 
to the quark mass difference rrid — m u , in order to make sure that isospin breaking 
effects do not produce unexpectedly large contributions. 

For a precision test of the unitarity of the CKM matrix, i.e. of the Standard 
Model, an accurate calculation of the RC, in particular a reliable determination of 
the effect of the hadronic structure, is crucial. The terms in the electromagnetic 
radiative corrections of 0(a) that are generated by the vector current (the first two 
terms in (1.2)) are firmly founded on a current algebra formulation and the details of 
the underlying quark structure are of only minor importance. We shall not further 
consider that part of the RC of 0(a) which is induced by the vector current. While 
the short-distance contribution of the axial vector current is well established too, 
its role at low energies strongly depends upon the detailed quark structure of the 
decaying hadron and its influence on the decay properties has been estimated only 
very roughly in terms of an effective low-energy cutoff Ma and the quantity C . We 
do not know of any published work that attempts to obtain the contribution of the 
axial vector current using a model of hadronic structure. However, it is evident that 
a reliable interpretation of the experimental data and a conclusive analysis of the 
unitarity problem necessarily requires a more refined treatment of the effect of the 
quark structure in order to substantially reduce the theoretical uncertainties and to 
firmly establish the size of the hadronic corrections. 

In this paper we shall calculate the axial vector contribution to the RC in the 
case of pion beta decay in the framework of the light-front quark model (LFQM), 
which is a relativistic constituent quark model based on the light-front formalism 
TBI. The LFQM provides a conceptually simple, phenomenological method for the 



determination of hadronic form factors and coupling constants, and has become a 
much used tool for investigating various electroweak properties of light and heavy 
(see e.g. [14, 15 1 and references therein). In Ref.|TJ| we have presented 



mesons 



a covariant extension of the LFQM which permits the calculation of all the form 
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factors that are necessary to represent the Lorentz structure of a hadronic matrix 
element. In this approach a meson is composed of valence quarks with constituent 
quark masses and the structure of the bound qq meson state is approximated by a 
covariant model vertex function, which depends on a parameter 1//3 which essentially 
determines the confinement scale, i.e. the size of the composite meson. Form factors 
are given in the one-loop approximation as light-front momentum integrals. As an 
example, it was shown in Ref. that a prediction of the electromagnetic form 
factor of the pion for small values of the momentum transfer can be made that is in 
good agreement with the data. 

The simple structure of the qq bound state should allow definite conclusions 
about the relative importance of the hadronic environment in a calculation of the 
RC. Radiative corrections of order a to the form factors that describe pion beta 
decay arise from the virtual exchange of Z, 7 or W and are represented by two- 
loop diagrams. We shall extend the approach of Ref. []16] and derive unique LFQM 
expressions (that are free of spurious contributions) for the two-loop diagrams asso- 
ciated with the axial vector current, and derive in this way the effect of the hadronic 
structure on the 0(a) corrections for pion beta decay. This determination of the 
effect of the hadronic environment by means of a two-loop calculation should be just 
as reliable as the one-loop calculation of the electromagnetic form factor of the pion. 

We shall show in this work that the uncertainty of the hadronic corrections due 
to the particular quark structure of the pion is small for pion beta decay. This 
result is in contrast to the situation for superallowed nuclear decays and neutron 
decay where the large value of C, Eq.(1.4), signals a much greater importance of the 
detailed quark structure with all its model dependent uncertainty. We shall analyze 
superallowed nuclear decays and neutron decay in a future work in a similar manner 
as for pion beta decay. But even without knowing the result of such an investigation 
it is clear that pion beta decay, once precise data are available, will always have 
a unique position due to the simple quark structure of the pion which generates 
hadronic corrections with very small uncertainties. 

In Sect. 2 we present the general formalism for pion beta decay without radiative 
corrections, which is analyzed in terms of two form factors that describe the quark 
structure of the pion. We investigate the effect on the decay rate of both the isospin 
violation due to the quark mass difference and the momentum transfer dependence 
of the form factors. In Sect. 3 the detailed calculation of the RC due to the axial 
vector current is presented. We consider three different groups of two-loop diagrams, 
and derive their light-front representations. We discuss the associated zero-mode 
problem in the Appendix and show that the respective light-front amplitudes are 
unique, i.e. free of spurious contributions. We approximate higher order gluon 
exchange effects by means of p exchange diagrams, which are shown to be of only 
minor importance if appropriate off-shell form factors are used. Sect. 4 contains our 
result for the RC for pion beta decay. 
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2 General formalism for pion beta decay without 
radiative corrections 

The amplitude without radiative corrections for the decay 7r + — > ii°ev is given 

by 

Ti = ^V ud ( n°(P")\d^(l - lM* + (P') )L» (2.1) 
where the matrix element of the leptonic current is 

Lfi, = u u {ku)l»{l ~ 7sK(0 (2.2) 

and k u , I are the 4-momenta of the neutrino and the positron respectively. We 
represent the hadronic matrix element for pion beta decay in terms of appropriate 
form factors 

(n°(P")\d lfl u\n + (P')} = ^{(P' + n.F^+q^q 2 )} (2.3) 

(n°(P")\d^ 5 u\n+(P')) = (2.4) 

where q = P' — P" is the 4-momentum transfer which varies within the range 
m 2 e < q 2 < (M+ - M ) 2 . 

It is convenient to analyze semileptonic decays of pseudoscalar mesons in terms 
of the form factors F\(q 2 ) and F (q 2 ), where the scalar form factor F (q 2 ) is defined 

by 

F °(' 2 > = F ^> + W^ F2(,2) ' (2 ' 5) 
The differential partial width in terms of these form factors is then 

rf r (7r+ -> 7i°eu) G 2 F \V ud \ 2 Ml 2 

dq~ 2 = 32.3 P(g ) ( 2 - 6 ) 

and the density p(q 2 ) consists of spin and spin 1 contributions as follows 

pfa 2 )=Po(« 2 )+P.fa 2 ) (2.7) 

^Hf(W))>-f)Vl) ! t <«> 

where p n is the recoil momentum of the n° in the ir + rest frame: 

p " = 4My{ (M + " M ° 2)2 + g4 " 2g2(M + + M ° 2) }- (2 - 10) 



If the quark structure of the pion is neglected, i.e. in the limit Fi(q 2 ) = F (q 2 ) = 1 
and for the approximation p 2 ~ (M+ + Mo) 2 (A 2 — q 2 )/4M+, the integrated partial 
width leads to the approximate expression for the total decay rate 1/ro, Eq.(1.13). 

In this Section we shall briefly discuss the exact integrated partial width (2.6) 
based upon the formulas for the form factors Fi(q 2 ) and F 2 (q 2 ), which we have 
derived in the framework of the quark model in Ref. | |16| |. 




P" 




(b) 



Figure 1: The one- loop contributions to pion beta decay. 



The hadronic matrix element (2.3) is given in the one-loop approximation, cor- 
responding to the diagrams of Fig.l, as a light-front momentum integral, denoted 
by A^. The 4-momentum of a meson of mass M' in terms of light-front compo- 
nents is P' = (P'~, P' + , P±), where the transverse vector is P' ± = (P n ,P' 2 ). Its 
constituent quarks have masses m^, m<i and 4- momenta p[, p 2 , respectively, and the 
total 4-momentum of the meson state is given by p[ + p 2 = P', i-e. the quarks are 
in general off the mass-shell. The appropriate variables for the internal motion of 
the constituents, (x,p' ± ), are defined by 



xP' + 
xP' ± + p' ± 



p+ = (l - x )p'+ 

, P2± = (1-X)P' ± - P ' ± 



and the kinematic invariant mass is 



p'l + m' 2 p'l + m\ 



x 



1 — X 



(2.11) 



For the transition between an initial 7r + = ud with 4-momentum P', mass M', 
and internal variables and masses of its constituent quarks (x,p' ± ,m' 1 ,m2) and a 
final 7T° = (dd — uu)j\[2 with 4-momentum P", mass M", and the corresponding 
internal quantities (x,p±,m",m2), the momentum integral A^, in a Lorentz frame 
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with q + = 0, consists of two parts that describe the u — > d transition of Fig. la and 
the d — > u transition of Fig. lb, and is given by 



A, 



1 

71 



H„(m„, m d , m d ) + H^(m d , m u , m u 



(2.12) 



where 



with 



H fJ ,(m l 1 ,m",m 2 ) 



16tt 3 



dx J d 2 p' L — 



h' h" 



x)N[N- 



Si, 



ii A 4 



= tr [75(^1 + m'Dj^^y + m' 1 )7 5 (- fa + m 2 )] , 



(2.13) 



(2.14) 

where N c is the number of colors, i. e. N c = 3. The light-front momentum integral 
H^, Eq.(2.13), is computed at the pole of the spectator quark: 



N 2 = v\ - m 2 2 = 0. 



(2.15) 



In our formalism [16[ four-momentum is conserved and the 4-vectors appearing in 
the trace (2.14) are then given by 



,m 



P2 



Pi 

P'~P2 



Pi 



(1- 



where m 2 _|_ = m 2 + p±. It follows from Eq.(2.16) that 



N[ = pf 



m' 2 = x(M' 2 - Mq ) 



P? 

J/2 



III 



n2 



112 



K 2 ) 



Ml 2 



P± + (1 — x)m" 2 + xmr 2 



x(l 



X 



(2.16) 



(2.17) 



and p'[_ = p' ± — (1 — x)q±. In our phenomeno logical approach we have chosen a 
pseudoscalar vertex operator for the qq pair bound in a S-state state, with the 



matrix structure of 75 and vertex functions h' and Ji'q, where [16 
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m 



2\2 



4M<f 



1/2 



M 



12 



M' 2 



[Mq 2 - (mi -m^ 2 ] 1 /: 



AO 



(2.18) 



for the qq bound state of mass M', and a similar equation for Hq. The orbital wave 
function is assumed to be a simple function of the kinematic invariant mass as 



0(Mq 2 ) = N'exp(-M' 2 /8/J 



ji2\ 



[2.19) 



where N' is the normalization constant and the parameter 1/(3' determines the 
confinement scale. The normalization condition is obtained for M' = M", (3' = 



f3" = (3, m; 



m 



m 2 = m and q 



0, either relation for H^(m,m,m): 



H^(m, m, m) 



{P 1 + P 1 ' 



;2.20) 
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or as a relation for the orbital wave function: 



which for the equal mass case is given explicitly by 

0( M( f ) = vr- 3 / 4 /T 3 / 2 (^) 1/2 exp( - (Mq - 4m 2 )/8/3 2 ). (2.22) 

While the form factor F\(q 2 ) in the one-loop approximation can be derived di- 
rectly from the plus component of the momentum integral (2.12), the calculation 
of the form factor F 2 {q 2 ) requires an appropriate account of the effect of zero-modes, 
as we have shown in []T6j . We shall not write down the formulas for the form factors, 
they can be found in Ref.[|l(|, but quote the results of the numerical calculation. In 
the limit of exact isospin symmetry the quark masses and the pion masses are equal, 
i.e. m u = rrid = m and M + = Mq = M n , and the form factors can be predicted: 
Fi(q 2 ) = F 1T {q 2 ), where F n is the charge form factor of the pion with i^(0) = 1, and 
F2{q 2 ) = 0. In our model the effect of isospin symmetry breaking is generated by 
a finite quark mass difference Am = — m u , while the parameters for the wave 
functions of 7r + and tc° are kept equal: f3 + = /3 = @ n . We use the parameters which 
we have found to reproduce the properties of pions in very good agreement with the 



data in Ref. IIB 



m = (m u + m d ) /2 = 260 MeV 

& = 308.8 MeV. (2.23) 

For the calculations of this Section we assumed a mass difference m^ — m u = 4 MeV. 

The momentum transfer in pion beta decay is small and the form factors can be 
approximated by monopole forms 



f ^ = Y^Ja 2 ' Z = ' L (2 - 24) 



The explicit calculation gives A x = 719 MeV (the corresponding quantity for the 
charge form factor of the pion is A n = 720 MeV), while 

Atyi 2 

Fi(0)wl- 7 - = 1 - 2.0 x 10" 5 , (2.25) 

v ; (902 MeV) 2 

from which it is seen that the effect of symmetry breaking on i*\(0) is of second 
order, in accordance with the Ademollo-Gatto theorem fI7 |. In contrast, ^(0) is of 



first order in the pion mass difference M + — Mq and takes the value 

F 2 (0) = -1.44 x 10" 3 , (2.26) 
which leads to the monopole approximation (2.21) for Fo(q 2 ) with Aq = 1.123 GeV. 
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The width Tq can be obtained from Eq.(2.6) by a numerical integration over q 2 
with the result 

r (7r+ -> 7r°ez/) = l/r (l - 1.2 x 1(T 5 ), (2.27) 

where 1/tq is the approximate expression given by Eqs.(1.13) and (1.14). The cor- 
rection is essentially due to the quark structure of the pion. Obviously, the effect 
of the symmetry breaking, Eq.(2.25), is largely compensated by the effect of the q 2 - 
dependence of the form factors, and the sum of all structure dependent contributions 
to the transition probabilities, Eq.(2.26), is indeed very small, and can be safely ne- 
glected. We shall continue to analyze pion beta decay in the isospin symmetry limit 
m u = = m, with m given by (2.23). 

3 The radiative corrections of 0(a) from the axial 
vector current 

The axial vector current essentially contributes to pion beta decay in 0(a) only 
in the two-loop processes which are represented by the vertex correction diagrams 
of Fig. 2 and the exchange diagrams of Fig.3. The amplitude corresponding to the 
photon-exchange diagrams of Fig. 2, involving the axial vector current, is given by 

7*> = %K. « fd<k ,„, .± L *l, , M "' . . (3.1) 



y/2 4tt 3 J (k 2 + ie)(k 2 - 2lk + ie) M 2 , - (q - k) 2 + is 
The leptonic tensor is 

Lfj,x = u u (k u )^^(l - 75) (-/+/; + m e )^ x v e {l) 

= -2/ A L M + A; A L M + A; M L A - g^kL + ie^k"^ (3.2) 

and the leptonic current has been defined in Eq.(2.2). The hadronic tensor 
contains only the axial vector part of the weak current. Current algebra methods 



have been used in Ref. |18| to derive its asymptotic behavior, which leads to the 
following expression for A^y. 

A, x = ~(Qu + Q d ) ie,x a pk a ( n°(P")\d^u\n + (P') ) - \ _ +0(-}-), (3.3) 



k 2 — M\ ^k 2 

where an arbitrary hadronic mass M4 is introduced to avoid a spurious infrared 
divergence in Eq.(3.1). The low-energy part of A^\ depends on the quark structure 
of the pion and is unknown. 

If the result (3.3) is inserted into Eq.(3.1) for and added to the corresponding 
Z-exchange contribution T 2 of Fig.2, one obtains the correction terms of 0(a) in 
Eq.(1.2) that are induced by the axial vector current, where the unknown low-energy 
contribution is parametrized in terms of the constant Ma- 

It is the main purpose of this work to evaluate those contributions of the vertex 
correction and exchange diagrams that come from the axial vector current, in the 



11 




Figure 2: Vertex corrections for pion beta decay. 



light-front quark model of Ref . [ffBl . The model calculation coincides with the result 
of Eq.(3.3) in the asymptotic limit and completes the current algebra approach by 
filling in the details that depend upon the quark structure of the pion. 

In addition we shall estimate the contribution of higher order gluon exchange by 
means of the p exchange diagrams of Fig. 4 

3.1 The vertex corrections of Fig. 2 

We shall calculate the contribution of the vertex correction diagrams of Fig. 2 
by treating separately the vertex correction for an off-shell quark. In the limit 
I = k v = q = the matrix element for the exchange of a photon, that consists of a 
part that describes the u — > d transitions of Figs.2a,b and an analogous part for the 
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d — > u transitions, is given by 

tf> = G , v- «=_ /; dx J ^j^—^ [gjiw + <M*»] (3.4) 

and the contribution of the vertex correction is 

n(°) = ^ [ d A k S $ — M ^ 

W (A; 2 + ie) 2 (A; 2 -2p / 1 A; + iVi + ie) M^-F + ie' 1 ; 

ttW = — />A; ^ L ^ M w (36) 

7T 2 J (k 2 + ie) 2 {k 2 + 2p'{k + Af + ie) - k 2 + *e ' 1 ' ' 

with 

tf$ = tr [75(^1 + K')(-7m75)(A- >* + m) 7 A(^i + ™ih 5 (- A + m 2 )] , 
S$ = tr [75 (A' + m'i)lx{£>'[+ A + ^)(-7m7 5 )(A + m / 1 ) 75 (- ^ + m 2 )] , 

where p'[ = p\ and N" = N[, since g = 0. The evaluation of the traces gives the 
result 



S$ = -ie^AN^pl + k a S?), 



where S@ has been defined in Eq.(2.14). Only the term ie^ Xa pk a L p of the lept onic 
tensor L^\, Eq.(3.2), contributes to the momentum integrals 11^ and IT^, which 
can be written as 

2 r . k 2 ■ LS - kS ■ kL + 4N[ (p[k ■ P'L - p[L ■ kP r ) 



n w = —(d'k 
in 2 J 



(k 2 + ie) 2 {k 2 - 2p[k + N[+ ie) M^-tf + ie 

(3.7) 



n w = —ld A k 

in 2 J 



k 2 -LS-kS-kL- AN[(p[k ■ P'L - p[L ■ kP') Mir 



(k 2 + ie) 2 (k 2 + 2p'(k + N" + ie) M 2 , - k 2 + ie' 

(3-8) 

where we have used that p 2 = P' — p[. The momentum integrals of Eq.(3.6) can be 
calculated in terms of the usual space-time components by the standard Feynman 
parameter method. Using the detailed results that have been collected in Appendix 
A we find 

n (o) = n (6) = n = 2{(3&! + b 2 p' 2 )LS - b 2 p[S ■ p[L 

-4ox N[{p' 2 ■ P'L - p[L ■ p^P')}, (3.9) 

where 01, b\ and b 2 are functions of p' 2 and are given by Eqs.(A3)-(A5). 
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The matrix element for the exchange of a photon has been expressed in terms 
of the light-front momentum integral (3.4) which is to be computed at the pole of 
the spectator quark. However, it is well known (see e.g. Ref.[|l6| and references 
therein) that this straightforward light-front representation of a hadronic matrix 
element is in general incomplete and contains spurious contributions that violate 
Lorentz covariance. These difficulties are a consequence of the fact that the effect 
of the associated zero-modes is not included. (Examples of zero-mode contributions 
can be found in Appendix B, Eq.(B5), and in Ref. ||16||. A more general discussion 
in the context of light-front quantization is given in Ref.[^9| .) This is the zero- 



mode problem which in the present case can be circumvented by the decomposition 
of the matrix element into a covariant (physical) part, that is not associated 
with a zero-mode, and a spurious part that is cancelled by the appropriate zero- 
mode contribution. We are only interested in the physical part of that can be 
identified by choosing a special representation of the 4-vector L: 

L = (L~, 0, 0_l). (3.10) 

The method which we have developed in Ref.[|l(| can be used to show that the 
resulting expression for T 2 (7) is unique, i.e. the contribution of the associated zero- 
mode vanishes exactly. 

We thus conclude that the condition (3.10) guarantees that all spurious contribu- 
tions are eliminated and the momentum integral (3.4), calculated at the pole of the 
spectator quark^ uniquely defines the complete light-front representation of the ma- 
trix element T 2 . In order to express the quantity II, Eq.(3.9), in terms of light-front 
variables, we compute the following scalar products, using Eqs.(2.15)-(2.17), 



p[L = xP'L 
4: 
1 



LS = 4xMq P'L, 



p[P' = ^M 2 n + N[), (3.11) 
p[S = 2(p?M 2 + m 2 M 2 -N[ 2 ). 

In this manner one finds II as a function of p' 2 = m 2 + N[ = m 2 + x[M 2 — Mq 2 ): 

IT = SP'L^hxMg - (ai + b 2 )N[ [pf - |(jVj + M 2 )] }. (3.12) 

Inserting (3.12) into Eq.(3.4) we can express the matrix element in terms of T\, 
Eq.(2.1), as 

T {l) _ l T Jf(27) 

1 2 — 2 1 axial ' 

- 2( Bl + b 2 )N[(p' 2 - X -{N[ + Ml))}. (3.13) 
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In arriving at Eq.(3.13) we have finally established the LFQM expression for the 
matrix element 

In order to complement the above remarks regarding the zero-mode problem, 
we note that the 3-dimensional light-front momentum integral (3.13) with pointlike 
nqq vertices (i.e. for h Q = const.), and the covariant 4-dimensional momentum in- 
tegral that represents the photon-exchange Feynman diagrams of Fig. 2 (with the 
same pointlike nqq vertices) are equal, which is another proof that there are no 
zero-mode contributions. We can reverse this argument and conclude that the tran- 
sition from the covariant Feynman perturbation theory to the LFQM proceeds in 
two steps: In the first step the manifestly covariant 4-dimensional momentum inte- 
gral, that corresponds to a given Feynman diagram, is represented exactly in terms 
of a 3-dimensional light-front momentum integral. In the second step appropriate 
phenomenological nqq vertex functions are introduced into the light-front represen- 
tation. 

For the amplitude of the Z-exchange diagrams of Fig. 2, that can be derived 
by an analogous analysis, one finds, that the quark structure contributes to the 
amplitude only to order m 2 /M^ and m 2 /M§; this is a small effect which can be 
safely neglected. Therefore, the published result, that is due originally to Sirlin [|TJ, 
remains essentially unchanged and is 

T (Z) _ 1 T A2Z) 
1 2 — 2 1 axial 

<£S = frQ^Hk + O^,*)), (3,4) 

where the Minimal Standard Model relation Myy = Mz cosQw (@iy is the weak 
angle) has been used (see e.g. Ref.@]). 

In the calculation of the vertex-loop integrals the value of the mass of the internal 
quark line is important only in the low photon momentum range, where the quark 
mass is essentially equal to the constituent mass. Therefore, for the numerical 
computation of Eq.(3.13), we take the values (2.22) for the parameters m and j3 n , 
and obtain the result 

A = -4.498, (3.15) 

where we have used the normalization condition (2.21). 
The combined correction is 

€m = *22, + «££ = tM q - + Q " )fn ii + \ + f )■ (3 - 16) 

3.2 The exchange corrections of Fig. 3 

The exchange diagrams of Fig. 3 are different from the diagrams of Figs.l and 
2 in that both quark lines are involved in the decay process, and there is no well 
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defined spectator quark line. In order to establish the LFQM expression for the 
amplitude T3 that corresponds to the exchange diagrams of Fig. 3 we shall use the 
analysis of the covariant two-loop diagram with pointlike nqq vertices ( a Feynman 
diagram) presented in Appendix B, guide. 



Just as in the case of the one-loop momentum integrals we can calculate the 
two-loop 4-momentum integrals, expressed in terms of light-front variables, by per- 
forming the integrals over the minus components of the loop momenta by contour 
methods, whereby the momentum integrals are given as residues of the respective 
quark poles. We have emphasized in the previous Subsection that this straight- 
forward procedure leads in general to an incomplete result since it misses the con- 
tribution of the zero-modes, and without including this effect the contour method 
not only violates Lorentz covariance, but is an uncertain approximation of the 4- 
momentum integrals. The zero-mode problem can be circumvented only if the spe- 
cial representation (3.10) is used. In Appendix B we prove that the contribution 
of the zero modes vanishes exactly for the resulting amplitude T 3 . Therefore, the 
contour method already leads to the complete result for the momentum integrals, 
which are given as residues of the respective quark poles, and uniquely defines their 
light-front representation. 

We shall calculate the amplitude T3 in the limit l = k u = 0, q=l + k ll = 
and P' = P", with P' 2 = M%. Only the contribution of the axial vector current 
will be considered. It consists of two parts, which depend upon the sign of the 
plus component of the photon momentum k = p'[ — p[ = p' 2 — p'l- For k + > the 
residue is determined by the poles of the quarks with momenta p' x and p'^. The 
remaining momenta are determined by 4-momentum conservation, i.e. p' 2 = P' — p[ 
and p'l = P' — p'l- These conditions lead to the relations: 




(a) 



(b) 



Figure 3: Exchange corrections for pion beta decay. 



k + > : 








n2 2 n 

P2 ~ m = °> 
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P'i-m 2 = (1 - x')(M 2 n - M?) , pf -m 2 = x"(M 2 -M{ 2 ), 

k 2 = kl = n{Ml - (1 - x')M^ 2 - x"K 2 } - hi, (3.17) 

where « = k + / P' + = x" — x' and < k < 1 — x'. 

For /c + < the residue is determined by the poles of the quarks with momenta 
p' 2 and p". These conditions lead to the relations: 

k + < : p'l - m 2 = , pf - m 2 = 0, 

pf-m 2 = x'(M 2 -M> 2 ) , pf -m 2 = (l-x")(M 2 -< 2 ), 

A; 2 = /c 2 = /t| — M 2 + x'M* + (1 - x")K 2 } - *a, ( 3 - 18 ) 

and /t varies within the range —x' < k < 0. The resulting amplitude T 3 is then 
given by 



h' h" 



x')x"(l - x")(M 2 - Mg)(M% - M<f ) 



(3.19) 

The hadronic tensors associated with the diagrams of Fig.3a,b are 

S$ = tr [ 75 (#' + m)(-7 M75 )(/ 1 + m) 75 (- + m) 7A (- ^ + m)] , 

45 = ^ [75(^2 + ^)7A(^2 + m) 7 5(- /i + m)(-7 M 75)(- it[ + m)\, (3.20) 



-5 



(a) 



The leptonic tensor L^x has been defined in (3.2). The product of the hadronic and 
leptonic tensors consists of two parts 

S^U* = ie^ al3 k a L p S { ^ + pseudoscalar, (3.21) 

where the pseudoscalar terms do not contribute to the amplitude T 3 and will be 
ignored in the following presentation. The scalar terms are given by 



S$L>* = S(p'[pl + m 2 ){p\k- P'L-p[L- P'k 

-8 (p[p 2 + m 2 ) {p'[k ■ P'L - p'[L ■ P'k^j . (3.22) 
Using the special representation (3.10) we find that 

S$V* = 8PiL(pi;p>> + m 2 )(p[k-x>P>k ) ) 

-8 P'L (p[p' 2 + m 2 ) {p'[k - x"P'k^j . (3.23) 
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This equation can be written such that its value at the various quark poles becomes 
obvious: 

S$U* = 2P'L(M 2 -(p'; 2 -m 2 )-(p'?-m 2 ))((l-x')(p'; 2 -m 2 ) 
_(1 _ x ')( p f _ m 2) + x yf _ m 2 } _ ^,2 _ m 2 } _ fc 2) 

+2 P'L (M 2 - (pf - m 2 ) - (p£ - m 2 )) ((1 - x")(p' 2 - m 2 ) 

-(1 - x"){j>'[ 2 - m 2 ) + x"(p' 2 - m 2 ) - x'\p" 2 - m 2 ) - k 2 ) . (3.24) 



The contribution of S$L^ to the integrand of (3.20) for k+ > is, according to 
the conditions (3.18), given by 

k + > : SftlJ* = 2 P'L- R(x', x"), (3.25) 

where 

R(x', x") = ((1 - x")M 2 + x"M'f) (x"(l - x')(M 2 - M'f) - x'{\ - x')(M 2 - M' 2 ) - k\ ) 

+ [x'Ml + (1 - x')M' 2 ) (x"{l - x')(M 2 - M' 2 ) - x"{\ - x"){Ml - < 2 ) - k\ ). 

(3.26) 

For k + < we use the conditions (3.19) to find 

k + < : SjjL"* = 2 P'L- R(l -x',1- x"). (3.27) 

The computation of the amplitude T 3 can be simplified by the observation that the 
integrals for k + > and k + < are equal, which can be shown by substituting x' 
for 1 — x' and for 1 — x". 

Inserting Eq.(3.26) into Eq.(3.20) we find for the LFQM expression for the matrix 
element T 3 

T.^-T^Li (3-28) 



where 

(Qu + Q d ) ^§ £ dx' J d 2 P ' L £ dx" J d 2 v \ 

h' hi R(x', x") 



x(3) 
axial 



x'(l- x')x"(l - x"){Ml - Mg)(M* - M^ 2 )ki 



(3.29) 



For the numerical calculation of 5^J ial we take again the values (2.22) for the param- 
eters m and f3 v and find the correction 

€L = £ 0.256. (3.30) 
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3.3 The p exchange corrections of Fig. 4 



Besides the diagrams of Figs. 2 and 3 there is also the sum of all irreducible higher 
order gluon exchange diagrams. The effect of this contribution can be approximated 
by means of diagrams of the type drawn in Fig. 4, where appropriate meson states 
are exchanged between the weak axial vector and the electromagnetic vertices. The 
Born approximation (exchange of a pion), which in general is expected to give the 
dominant contribution, vanishes, since the pion does not couple to the axial vector 
current. We shall consider only the lowest mass exchange process, i.e. the p exchange 
diagrams of Fig.4, and calculate the corresponding amplitude first for on-shell vertex 
structures, which are defined by the appropriate matrix elements. 





(a) 



(b) 



Figure 4: p exchange corrections for pion beta decay. 



The matrix element of the electromagnetic current for the transition p° — > 7r° is 

h = QuU^xu + Qddjxd 



(P'\j x \P";U 3 ) = eT 



(7) 
A;/ 



= e v {Qu + Qd)g{k 2 )ie XvaP P a k p . (3.31) 
The matrix element of the axial vector current for the transition 7r + — > p° is 



(P"; 1J 3 | - d W u\P') = >/2e*T 



*i/-p(A) 



= V2e*» { - f(k 2 )g^ - a + (k 2 )P^ + a^k 2 )^}, 

(3.32) 

where e = e( J3) is the polarization vector of the p, P = P' + P", and k = P" — P'. 

The amplitude that corresponds to the diagram of Fig.4a is then given in the 
limit I = k v = by 



T^a — —G V u( i{Q u + Qd) 



d A k- 



p(T) n upy(A) j u\ 

\v y up 



(2vr) 4 J ™ '" (k 2 + te) 2 ((P' + k) 2 - M 2 + ie) 
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(3.33) 



where we have used the bare p propagator 



Ao 1 ^) 



9»u 



P P 

P 2 



MP), 



P 2 - M 2 p + ie, 



(3.34) 



and the leptonic tensor is given by Eq.(3.2) and P' 2 = M 2 . Using the hadronic 
tensors as defined by (3.31) and (3.32) gives the result 



T^a — —GV u d{Q u + Qd) 



ie 



(27T) 



d±kg{k 2 )f{k 2 



4(P'k -kL-k 2 ■ P'L) 
{k 2 + ief(JP' + k) 2 - M 2 p + ie) ' 

(3.35) 

In the isospin symmetry limit the amplitudes corresponding to the diagrams of 
Fig. 4a and 4b are equal, and the total contribution is given by 



Ta — Td a + T z 



-46 



2T Z 



4a • 



The amplitude T4 depends only on the form factors g(k 2 ) and f(k 2 ), which we have 
determined in the framework of the light-front formalism in Ref . |16| . The results of 



16 can be written as 



K{M' Q )h p {MH) 



x < m + 



- x)x 2 (P> 2 - M<?)(P" 2 - Mq 2 ) 

(p' ± k ± r 



Mq + 2m 



V'l 



k 2 

K{M' Q )h p {M{) 



(3.36) 



Ii2\ 



X 



- 2xmM'£ ■ 



a;)a; 2 (P /2 -M^)(P" 2 -M ( 



2xmP' 2 - mkP + mk 2 + 2m(£; 2 - kP) 



Mq + 2m 



2xP' 2 + 2xM' Q 2 — k + kP — 2{k - kP) 



P±k± 
k 2 

V'\_k\ 



(3.37) 



where p'[_ = p'j_ — (1 — x)k_L, k\ = —k 2 and kP = P" 2 — P' 2 . We have designated the 
nqq vertex function h' by h^{M'^) and the pqq vertex function by h p (M{{); they 
are given by Eq.(2.18) in terms of the 7r mass M v and p mass M p , respectively. The 
on-shell form factors are given by 



g(k 2 )=g(k 2 ;M 2 n ,M 2 p ) 



f{k 2 ) = f{k 2 -MlM\ 



(3.38) 



For the evaluation of the 4-momentum integral (3.35) it is convenient to approximate 
the form factors by monopole forms: 



g(k 2 



9(0) 

-k 2 /A 2 



f(k 2 



IB 

-k 2 /AV 



(3.39) 
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where the pole masses A g and A/ are determined by the derivatives of the form 
factors at k 2 = 0. For the numerical calculation we take the values (2.22) for the 



parameters m and (3^, and f3 p = 0.26 GeV [flql , and obtain the results 

3(0) = -1.21 GeV' 1 , A g = 0.6G4GeV, 

/(0) = -0.85 GeV , A f = 1.72 GeV. (3.40) 

The 4-momentum integral (3.35) can now be calculated by the standard Feynman 
parameter method, with the result 

T - l T A (4) 

SaL = |-(-0.69) = -8.0 x HT 4 (on-shell). (3.41) 

Sirlin has estimated the contribution of the diagrams of Fig.4 on the basis of vector 
dominance and Weinberg sum rule arguments in Ref. |J and found the correction to 
the decay rate 'to be a few times 10~ 4 ', in accordance with (3.41). 

However, the intermediate p in the diagrams of Fig.4 is off-shell, and for a rigor- 
ous evaluation of the corresponding amplitude the off-shell structure of the hadronic 
vertices must be accounted for. Moreover, a consistent treatment requires the cal- 
culation of the p self-energy operator in the same gg-loop approximation that has 
been used for the calculation of the hadronic form factors. The corresponding renor- 
malized p propagator is obtained from (3.34) by the modification 

Ao \P) - A- X (P) = (P 2 - M 2 + is) P p (0; M 2 , P 2 ), 

where F p (k 2 ; M 2 , P 2 ) is the half-off-shell charge form factor of the p, with the nor- 
malization condition P p (0;M 2 ,M 2 ) = 1. An analogous result has been derived for 
the renormalized pion propagator in Ref.[[2(|. For k 2 = the charge form factors of 
7r and p are given in the one-loop approximation by the same analytical expression 
(with appropriate 7r and p parameters) which has been derived in Ref. [T^] to be 



(3.42) 



In order to estimate the importance of this structure effect we have used the light- 
front formulas (3.36) and (3.37) to continue the form factors g(k 2 ) and f{k 2 ) to their 
off-shell forms g(k 2 ; M 2 , (P' + k) 2 ) and f(k 2 ; M 2 , (P' + k) 2 ), with P' 2 = M 2 . For 
an order of magnitude estimate it is sufficient to compare the values of the different 
form factors at k = 0. We find 

g(0;M 2 ,M 2 ) = 0.125 GeV' 1 , /(0;M 2 ,M 2 ) = -0.028 GeV, P p (0;M 2 ,M 2 ) = -0.070, 
and consequently 

g(0; M 2 , M 2 )/(0; M 2 , M 2 )/(^(0)/(0)P p (0; M 2 , M 2 )) = 0.048. (3.43) 
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Therefore, without going into any further computational details one can conclude 
that the resulting correction is reduced by at least the factor (3.43), i.e. 

Cai = O(l0- 5 ) (off-shell). (3.44) 

Evidently the correction due the p exchange diagrams of Fig.4 are very small, and 
can be safely neglected. 



4 Concluding remarks 

We have developed in this work a strategy to deal with the effect of the hadronic 
structure on the corrections to pion beta decay in 0(a) due to the axial vector 
current. It is based on the light-front quark model for the pion, whose qq bound 
state structure is well described by two adjustable parameters, constituent quark 
mass and confinement scale, that have been shown in Refs.||14||, [ |I~6|j to describe a 



large body of data. In most applications of the light-front quark model the effect of 
the hadronic structure could be well approximated by one-loop diagrams. However, 
the axial vector 0(a) corrections to the beta decay of a qq bound state requires the 
calculation of two-loop diagrams, and we have considered three different types which 
are represented by the graphs of Figs. 2, 3 and 4. The corresponding amplitudes can 
be expressed in terms of light-front momentum integrals, which we have shown to 
be unique, i.e. there are no associated zero-mode contributions. 

The results of the quark model calculation for the corrections of 0(a) due to 
the axial vector current are given by Eqs.(3.16), (3.30) and (3.44), which can be 
rewritten in the form of the standard representation (1.2) in terms of the quantities 
M A and C: 

€L + €L = ^3(Qu + Q d )in^, 
Z7T Ma 

M A = 425.68 ±8 Mel/, (4.1) 

and 

*£L = ^(Qu + Q d )2C, 

C = 0(l(T 2 ). (4.2) 

Comparison with (1.3) shows that the value (4.1) for the effective mass M A is clearly 
different from the a\ meson mass, but close to the confinement scale of the qq pion, 
and is even below the guessed range (1.3). The error bar associated with Ma in 
Eq.(4.1) is due to the small uncertainties of the constituent quark mass m = 260 ±5 
MeV and the Z mass M z = 91.188 ± 0.007 GeV §. The value of C, i.e. 
the correction due to the p exchange diagrams of Fig.4 is so small that it can be 
neglected. The same is true of the effect of the form factors of the pion which we 
have discussed in Sect. 2. 
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Thus our model calculation of the axial vector contribution to the radiative 
corrections of 0(a) not only gives a definite value for Ma, Eq.(4.1), which leads 
to larger corrections, but also removes the large uncertainty in the RC due to the 
assumed range (1.3). Using the average value of Ref.|| A = M + — M = 4.5936 ± 
0.0005 MeV, which gives the end-point energy E = (M + + M ) A/2M+ = 4.5180 ± 
0.0005 MeV, we find (a/2n)g(E ) = 1.0515 x 10" 2 and obtain from Eqs.(1.2) and 
(4.1) the value of the RC to pion beta decay 

5 = (3.230 ± 0.002) x 10" 2 , (4.3) 

where the error comes from the uncertainty of Ma, Eq.(4.1). We emphasize that 
the error on the RC is of the same order as the neglected corrections (2.27), due to 
the weak form factors of the pion, and S^j ial , Eq.(3.44). 

We shall investigate the effect of the detailed structure of the three quark bound 
state on the properties of superallowed nuclear decays and neutron decay in a future 
work, in a similar manner as for the pion. If we tentatively assume that the effective 
mass Ma associated with a qqq nucleon will be found approximately equal to the 
result given by Eq.(4.1), it is easy to show that the unitarity sum derived from 
nuclear decays becomes V 2 = 0.9956 ± 0.0011, i.e. the violation of unitarity seems 
to be much more pronounced than indicated by Eq.(1.6). 

However, a definite judgement of the unitarity problem will be possible only if 
the complete results of both the detailed structure calculation for qqq nucleons and 
the precision measurement of pion beta decay will be available. 
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A Appendix: The momentum integrals for Sect. 3.1 

The momentum integrals encountered in the calculation of the vertex correction 
diagrams of Fig. 2 can be evaluated by the standard Feynman parameter method. 



Convenient formulas for the on-shell case can be found e.g. in Ref.[^l|]. In Sect. 3.1 
the vertex correction is calculated for off-shell quarks in terms of the function II, 
and the relevant integrals are given below: 



1 r ... k„ M 2 



-•-2 / ^ ( k 2 + i£ y( k 2 ± 2p[k + N[ + is) M 2 v -k 2 + IE ±P '^ aU (A ' X) 



iir 2 



d% {k 2 + ie) 2 {k 2 ±% i k + N[+ie) Nfy-V + ie ~ 9i * M +P >'^ (A ' 2) 
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where 

1 m 2 „ m 2 — p'? „ / 1 \ ... 

4 m J 8 4p'f 4p x 4 m 2 V M^ / / 

62 " 2pf p$* m pf iU m 2 (A ' 5) 

B Appendix: The two-loop diagram and its zero- 
mode contribution 

We shall analyze in this Appendix the two-loop diagrams of Fig. 3 for the special 
case of pointlike nqq vertex functions. From the corresponding covariant amplitude, 
which is given by the Feynman rules, we shall derive the light-front amplitude by 
integrating over the minus components of the momentum variables. For the purpose 
of this calculation we put the 7rqq coupling constant equal to 1. In the conventional 
space-time formalism the covariant amplitude in the limit I = k u = is given by 

T Feynman _ 1 n T/ ^^^c f ,4 / f ,4 // (g^M ~ Qi^jiS)-^ 

13 ' ~ 2 (2vr)8 J d Pl J d Pl D' x D' 2 D'lD'i (k 2 + is) 2 

X M*,-k* + ie' { ] 

where D' n = p' 2 —m 2 + ie, D" n = p'^ — m 2 + ie for n = 1,2, andp^+p^ = Px+p'i — P' 
with P /2 = M. 2 . The photon has momentum /c = p'[ — p[ = p' 2 — p 2 - The leptonic 
tensor has been defined in (3.2), and the hadronic tensors S$ and in (3.20). 
If the momenta are decomposed into their light-front components we have 

d 4 P [ = ^p' + dp'rdx'd 2 p' ± , d* P 'i = l -p !+ dp'[-dx"d 2 p" ± . 

We are only interested in the integration with respect to p'{~ and p'[~ , and use the 



same technique as in Ref.|]l6j, which is based upon the integral representation 

= f°° dae ia{p2 - m2+ie \ (B.2) 

p2 _ m 2 _|_ ^ J Q 

A similar procedure has been used in Refs.|22|] to investigate the relation between 
the standard covariant quantum field theory and light-front field theory. 

There are three basic integrals that contribute to the amplitude (Bl); these are 



dp'r J dp'; 



1 M 2 , 



2tt J J J D' X D' 2 D'(D' 2 ' (k 2 + is) 2 - k 2 + xe 
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Ml x'(l - x')x"(l - x")(Ml - Mg)(Mt - Mq 2 ) 
A Q{x"-x>) Ml Q(x'-x") M^ w } 

i ^ mi-% + k% Mi-k\V [a - 6) 



% V l A I" l A Pi Ml 

' Pl J Pl D' X D' 2 D'[D' 2 ' (k 2 + ief Ml-k 2 + is 



1 



Ml x'(l - x')x"(l - x"){Ml - M(?)(Ml - M 



f (Ml-x>M> 2 )Q(x"-x>) Ml {l-x')M' 2 Q{x'-x") Ml ) 

I ki Ml-ki + ^ Ml-k\Y K * A) 

where k\ and k\ have been defined in (3.17) and (3.18), and we have used that 
P' + = M n . The result for the first two integrals (B3) and (B4) coincides with the 
result obtained by the contour method, i.e. in both cases one finds the residues of 
the respective quark poles, and zero-modes do not contribute. 
The third integral we represent as 

-i-V fdv" fdv"- pTp ' r M * w - 

2ti J J Pl J Pl D'^D'lD'i (k 2 + is) 2 Ml - k 2 + is ~ 

-^{PW ~ A + P^{x' - x") + R 1 5(x')S(x") + R 2 S(1 - x')6(l - z")}- (B.5) 

The form of the residue terms is obvious 

(Ml-x'M' 2 ){l-x")M^ 



Pi = 



x'(l - x')x"{\ - x"){Ml - M' Q 2 )(Ml - Mq 2 ) ki ' 



p {l-x')M' 2 {Ml-x"M^) 
2 ~ x'{\ - x')x"{\ - x"){Ml - Mq)(MI - Mq 2 ) k^ { ' ' 

The functions R n = Rn(p' ± ,p±) for n = 1,2, which determine the zero-mode con- 
tribution, are independent of M n , but their detailed form cannot be derived by the 
method used above. For example, Ri is found to be the ratio of two functions, each 
of which becomes zero for x' — x" — 0. However, there is an alternative way to 
determine R\ and R 2 - It uses the limiting behavior of the integral 



ton, / <i% J m ^f— ? = O(MJ), (B.7) 



which is implied by Lorentz covariance. 
If we define 



P, (0) = lim P 1 = -, 

1 Af„-o (l-x')x"ki' 

P, (0) = lim P 2 = — 1 -— , (B.8) 

2 M^O 2 X'{l-X")ki' y ' 
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and integrate (B5) with respect to x' and x", then, according to (B8) the contribution 
of 0(M~ 4 ) must vanish exactly, which gives the conditions 



f 1 dx' f 1 dx"Pi 0) + R 1 = 0, 

JO Jx' 

L 



dx' / dx"P 2 w + R 2 = 0. (B.9) 



o 

From the Eqs.(B9) for i?i and R 2 it can be seen that R\ = R 2 , by substituting x' 
for 1 — x' and x" for 1 — x". 

This derivation of the zero-mode contribution shows clearly that a residue term 
which is derived by the contour method, may contain a spurious part that is not 
consistent with the requirement of Lorentz covariance. It is the zero-mode contribu- 
tion that cancels this unphysical part of the residue term. This is an example of the 
deep connection between the zero-mode and the Lorentz invariance of the light-front 
formalism. 

Next, we shall decompose the product of leptonic and hadronic tensors in the 
integrand of (Bl) into products of light-front components. We use that = — S 1 ^, 
Eq.(3.20), and the result (3.23), and find 

S$L^ = 2P'l{ - 2 P '-p'[-{k+f 

+ p'f [k + {2m"l + x"Ml) - (1 - 2x')M n (2p'[k ± + x"M n k + 
~k + (2m'l + x'Ml) - (1 - 2x")M 7T (2p' ± k ± + x'M w k + ) 



- V"{ 



+ . 
(B.10) 



where we have omitted all those terms that are independent of p'i and p'{~ . 

If the decomposition (BIO) and the basic integrals (B3)-(B5) are used to perform 
the integration of (Bl) with respect to p'{~ and p"~, it is obvious that the zero-mode 
contribution of Eq.(B5) vanishes, since the term that contains the product p'{~Pi~ 
is multiplied with k + . The result for the momentum integral (Bl) is given by the 
residues of the respective quark poles, and zero-modes do not contribute: 

rpFeynman = Q ^ p , L ^ + ^ ^ ^ J ^ £ ^ J ^ 

R(x', x") 

x'(l - x')x"(l - x")(Ml - M' Q 2 ){Ml - Mi; 2 ) k% 
M 2 

(BU) 

where we have used that the integrals for k + > and k + < are equal, and the 
function R(x',x") has been defined in Eq.(3.26). 

We have proven that the covariant Feynman integral (Bl) and the light-front 
integral (Bll) are equal. Since the intgrals are finite, this result can be verified by 
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numerical calculation. At this stage we depart from the covariant (Feynman) per- 
turbation theory and introduce phenomenological vertex functions, Eq.(2.18), into 
the two-loop light-front integrals. This step gives the light-front quark model expres- 
sions, Eqs.(3.28) and (3.29), for the amplitude T 3 , corresponding to the diagrams of 
Fig. 3, in terms of a simple convolution of light-front vertex functions. 
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